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Abstract 



We investigate the heavy baryons which arise as solitonic excitations in a 
"heavy meson" chiral Lagrangian which includes the light vector particles. 
It is found that the effect of the light vectors may be substantial. We also 
present a simple derivation which clearly shows the connection to the Callan- 
Klebanov approach. 
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1 Introduction 



The problem of describing baryons containing heavy quarks (c, b, ...) as 
Skyrme type solitons may be considered as an extension of the problem of 
describing strange "light" baryons. An approximation in which the strange 
quark is considered very heavy was actually discussed some time ago by 
Callan and Klebanov 1 and others 2 . In this picture one imagines a sta- 
tionary K-meson to which is bound a nucleon treated as a soliton of the 
SU(2) x SU(2) effective chiral Lagrangian. Of course it is not clear that 
the approximation of treating the strange quark as "heavy" is the best one; 
comparable results may be obtained in other treatments 3,4 ' 5 which consider 
all the baryon octet members on the same footing. On the other hand, it is 
reasonable to expect that the bound state approach is a suitable one for the 
baryons containing c and b quarks. 

An application to the heavy quark baryons based on a fairly literal ex- 
tension of refs. 1 — 2 has in fact been given 6 . More recently it has been 
realized that to go to the heavy quark limit one should do more than just 
let the heavy mass get large. One should 7 at the same time keep the heavy 
particle four velocity fixed and impose the Isgur-Wise heavy spin symme- 
try. This has the consequence that a heavy vector particle should always 
appear together in a multiplet with the heavy pseudoscalar. The latter fea- 
ture was not taken into account in refs.l — 2, 6 and so the problem was 
recently reexamined in some papers by a San Diego-Caltech group 8,9 and by 
a Hanyang-Seoul- Jagellonian-Saclay-Stony Brook group 10 ' 11 . 

In the present note we shall consider the problem when the light vectors 
are included in the chiral Lagrangian. It is known that this makes both the 
mesonic and solitonic sector more realistic for the light particles. It is also 
likely to be important for the interaction of the heavy particles with the light 
ones; for example, the semileptonic D — > K* transition appears to dominate 
over D — > Kit. 

In addition, we will give a method of calculation which appears a bit 
simpler than the previous ones and which makes clearer the connection with 
ref. 1. This may be of relevance since there is apparently some disagreement 
between the two groups mentioned above. 
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2 Effective Lagrangian 



The light part of the effective Lagrangian of pseudoscalars and vectors, Cu g ht 
which we shall employ has been discussed in a number of places; see ref. 12 
for details and ref. 13 for an updating of the symmetry breaking part as well 
as further references. The relevant light fields belong to the 3x3 matrix of 
pseudoscalars, and to the 3x3 matrix of vectors, p M . It is convenient to 
define objects which transform simply under the chiral group: 

e = exp(^/F w ), U = i 2 

F^{p) = d^pv - d v p^ - ig [p„, p„] , (2.1) 

where F n m 0.132 GeV and g ~ 3.93 for a typical fit. We can do without 
the explicit form of Cu g ht here; it is given in (2.6), (2.7) and (2.13) of ref.12 
(wherein g should be replaced by g). We also do not need the details of 
the SU(3) symmetry breaking discussed in ref. 13. Actually we will restrict 
ourselves to two flavors of light quarks for simplicity. 

The new feature associated with the interactions of the heavy meson fields 
involves the use of a field 14 combining the heavy pseudoscalar, P' and heavy 
vector, Q' at fixed 4-velocity, : 

H= l - l ^ V » {i lb P' + i lv Q' v ), # = 74^74- (2.2) 

Note that H is here considered to have dimension of mass. H is a 4 x 4 
matrix in the Dirac space and it also carries an unwritten light flavor index. 
The chiral interactions of H with the light pseudoscalars were discussed in 
ref. 15. The extension to including light vectors was given in ref. 16, which 
notation we precisely follow here, and in ref. 17. The total light-heavy 
interaction terms we shall require are 16 : 

h ™ VV = iVf.Tr \H(d IJi - iagp^ - i{\ - a)v^)H] + idTr \H^^ b p^H 



%C r — "i 

+ _ T r \H lfllu F^(p)H] , (2.3) 
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where m v ~ 0.77 GeV is the light vector mass and 

v„p,= Z -(^±^d^)- (2-4) 

Furthermore M is the heavy meson mass, a, c and d are dimension- 
less coupling constants for the light-heavy interactions. The choice a = 1 
corresponds to a natural notion of light vector meson dominance, which is 
interesting to test. As seems appropriate for an initial treatment, we do not 



include terms in ( |2.3| ) which are higher order in 1/M, contain more deriva- 
tives of the light fields or involve light flavor symmetry breaking. 

For treating the light baryon as a soliton which gets bound to the heavy 
meson, we shall need information about the classical soliton solutions of 
Cught- First note the 2x2 matrix decomposition of p^. 

p, = + r>;). (2.5) 

Then we have (see e.g. ref. 12) the classical "profiles" : 

U c = exp(ix ■ rF(r)), u) 0c = Lj(r), 

Pic = ■^r^ika^kGir), (2.6) 

and uj ic = Pq c = 0. The boundary conditions for a finite energy baryon 
number one solution are: 

F(0) = -7t, G(0) = 2, u/(0) = 0, 

F(oo) = G(oo) = cu(oo) = 0. (2.7) 

Typically, the opposite sign is taken for F(0) as a kind of convention but 
the sign choice above gives the correct sign for both axial vector current and 
vector current matrix elements 18 . 



3 Overview of the approach 

The first step in the bound state approach is to find the channel in which the 
potential between the heavy meson and the nucleon as soliton is maximally 
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attractive. Of course in order to get a positive parity heavy baryon we should 
consider the orbital angular momentum I to be unity. In ref. 1, the attraction 
is found in the state with / = 1 and grand spin, G = |. The grand spin is the 
sum of the isotopic spin of the heavy meson, I and the angular momentum, 
J; it enters in the first place because the Skyrme soliton ansatz, given in 



(|2.6| ), mixes the isospin and ordinary spaces. We would like to display a 



suitable wavefunction for the heavy meson, H in the I = 1, G = | channel. 
First note from ( |2.2j ) that when we specialize the heavy meson to be at rest 
(Vi = 0), the 4x4 matrix H has non- vanishing elements only in the lower 
left 2x2 subblock. The first index of this submatrix represents the spin of 
the light degrees of freedom within the heavy meson while the second index 
represents the spin of the heavy quark. Furthermore there is the unwritten 
bivalent isospin index. Specializing to this subblock we then write the H 
wavefunction as: 

5* =(jfc. T ) M -^=-^L (3.1) 

Here, the radial wave function, u(r) is taken in the first approximation to be 
localized at the origin, r 2 u 2 (r) rs S(r). Note that the quantity x represents 
the angular part of the space wave function and the first factor couples it to 
the isospin index, b to give G = |. In turn this is coupled to the light spin 
index, I with the Clebsch-Gordan coefficient -^id to give G = 0. Finally, 
the heavy spin index, h is left uncoupled (as appropriate to the heavy spin 
symmetry) to give the desired net result G = |. 

To see that the wavefunction (|3.1|) gives attraction we consider the matrix 



element of the "potential" obtained from ( p.3| ): 

V = -iMd J d z xTr(H lilbPi H) + ■■■, (3.2) 

where only the d term has been shown for simplicity. Substituting in ( |3.1| ), 
Pi from (|2.4|) and (|2.6| ) as well as 7^75 — * — ioi yields 

V = ^F'(0)e sc (ai) s i j(x • T) ca [-^T; + x { ± ■ T] ab (x • r) bd j e u H . (3.3) 

Noticing that the object within the curly brackets is simply \{Ti)cd we easily 
get 

V = -yF'(O) + • • , (3.4) 
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This provides an attractive force in the I = 1, G = ~ channel since F'(0) > 
and d is also expected to be positive. We will give the contributions to V 
from the first and third terms in ( |2.3| ) in the next section; they do not appear 
to change the attraction to repulsion. However if we were to replace the 
Clebsch-Gordan coefficient in (|3.1| ) by SuSai, which corresponds to coupling 
everything but the heavy quark spin to G=l, we would find that the — | 
factor in ( p.4|) is replaced by +| and repulsion. This is in agreement with 
refs. 8 — 9; in our approach the x- r factor appears in a very natural way and 
also there is no need to introduce a collective rotation factor at the present 
stage. 

The next step is to introduce collective coordinates by allowing those 
constant transformations which leave the Lagrangian invariant to depend on 
time. First consider the soliton "rotational" modes, 19 A(t); these are defined 
from the transformations 

C/(x,t) = AtyUcWA- 1 ®; t ■ Pi (pc,t) = A(t)r ■ p ic (x)A-\t); 

H(x,t)=A(t)H boundl (3.5) 

where H bound is effectively the wavefunction in Q3.1]) and A(t) acts on its 
isospace index. Defining the angular velocity, f2 by A^A = |f2 • r and 
substituting (3.5) into / d 3 x(£u g ht + C-heavy) gives an additional contribution 
to the Lagrangian of the general form 1 : 

5L = —K Q 2 — x & ■ G H , (3.6) 

in which K and \ represent spatial integrals over the profiles in ( j2.6| ) while 
G H , to the accuracy of interest, is the expectation value of the heavy particle 
grand-spin, G H = (r/2) + • • • . Note that at the collective Lagrangian level, 
as discussed in ref. 1, the angular momentum of the heavy field is represented 
by its grand-spin while the heavy field isospin is transmuted to zero. From 
(|3.6|) , with the soliton angular momentum, J? = d5L/dQi, the rotational 
collective Hamiltonian is obtained as 

^4( JS + ^ fl ) 2 ' ( 3 - 7 ) 

The moment of inertia, K is identified from the light soliton sector as K = 
§(n"i A — ^iv) -1 in terms of the nucleon and A masses. Eqn ( |3.7| ) can be 



6 



simplified by noting that the total angular momentum is given as J = J s + 
G H . Then we arrive at the heavy baryon mass formula 1 : 

H coU = i(m A - m N ) [(1 - xW + 1) + X J(J + 1) + • • •] , (3.8) 

where the three dots stand for a term which does not split the heavy baryon 
masses. In deriving ( |3.8j ), use was made of the fact that the soliton carries 
the full isospin of the heavy baryon at the collective level so that J 2 = I 2 . We 
will see that Cheavy in ( P-3| ) leads to x = 0. Thus we have the final results: 

m(Z* Q ) - m(E Q ) = 0; m(Z Q ) - m(A Q ) = jj(m A - M N ) , (3.9) 

wherein the subscript Q denotes the heavy baryon in which the s quark has 
been replaced by the heavy quark Q. The equality m(Eg) = ?tj,(Sq) is, of 
course, expected from the heavy quark spin symmetry. 

Finally, we take into account the translational mode, R(t) for the soliton 
by now setting 

f/(x, t) = A(t)U c (x - R(t)) A-\t) , (3.10) 

and similarly treating the light vector fields. Imagining the heavy meson at 
rest with the soliton moving about it, we have a Schrodinger equation for the 
relative motion 

(~2^b v * + vm ) u(R) = Eu(R) ' (3,11) 

wherein V(R) is to be identified with ( |3.4|) . It will be expanded to quadratic 
order 

V(R) = V + ^kR 2 --, (3.12) 
following the notation of ref. 9. 



4 Effects of light vector fields 

Here we give details of the calculation, including all three terms in Cheavy, 



( |2.3|) . The a term involves the profile ui{r) and the c term involves the profile 
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G(r). The needed factors for the non zero contributions to the "potential", 
V(r) are 



1 

V2 



u(r), pi 



sin F 
2r ' 



1 



Xi x • r 



F'x,- x ■ r, 



-G'tr ■ t + ~{rG' + G(G - 2))x • a x • r 



We expand the profiles around the origin as follows: 

F(r) = -7T + rF'(O) + -r 3 F"'(0) + 

6 



(4.1) 



G(r) 



^r) = w(0) + -rV'(O) + ■ ■ 
= 2 + W(0) + ^r 4 G""(0) + 



(4.2) 



where powers were deleted so that the underlying fields are analytic 9 in the 
Cartesian coordinates. Substituting ( |4.1| ) and (JO) into ( |2.3| ) and evaluating 
the resulting expression in the attractive channel corresponding to wavefunc- 
tion (|3.1|) enables us to identify the parameters of the "potential" as 



V = ~dF'(0) + *^G"(0) - ^L(O), 
2 m,,^ V2 



(4.3) 



k = d(- 5 -F"'(0) + iF'(O) 3 ) + ^(|g""(0) + iG"(0) 2 ) - ^L"(0). (4.4) 
o 3 m„(7 o 2 v2 

If we choose to evaluate V(r) in the "repulsive" channel mentioned after 
( |3.4|) , the c and d terms in both ( |4.3D and Q4.4Q would be multiplied by —1/3 
while the a terms would remain unchanged. 

It is of course necessary to know the derivatives of the F(r), G(r) and 
cu(r) profiles which appear above. These are found by solving the coupled 
differential equations which arise from minimizing the static energy of £u g ht 
(see (6.6) of ref 12) subject to the boundary conditions ( |2.7| ) above. For 
the typical best fit to meson and light baryon masses these are uniquely 
characterized by the starting values 20 : 



F'(0) = 0.795 GeV, G"(0) = -0.390 GeV 2 , w(0) = -0.094 GeV. (4.5) 
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From ( |4.5| ) and the differential equations it is straightforward to deduce the 
next higher derivatives: 

F"'{0) = -0.12GeV 3 , G""{0) = 0.35 GeV\ u"{0) = 0.016 GeV 3 . (4.6) 

Finally, let us justify the vanishing of the \ f2 ■ (r) /2 term in (|3.6|) which 
was used to obtain the mass formulas in ( |3.9| ). This term might arise when 
( |3.5|) is substituted into (|2.3| ). We note that a possible contribution from 
the d term would involve a factor and hence a factor 7475, which has no 
components in the 2x2 subblock for which HH is non-vanishing. Similarly 
the relevant factor for the c term involves a factor 7^74, which also has a 
vanishing overlap with HH. For the first term in ( |2.3| ) we may explicitly 
check that the operator r has zero expectation value in the wavefunction 

5 Results and discussion 

To gauge the accuracy of the heavy quark symmetry we may observe that 
the vector-pseudoscalar mass differences,which should vanish in the M — > 00 
limit, are respectively about 400 MeV, 145 MeV, and 45 MeV for K* - 
K,D* — D and B* — B. From this point of view, the s-quark should not be 
considered heavy while it seems quite reasonable to treat the 6-quark in this 
manner. 

First let us look at the mass splitting relations (3.9), which are indepen- 
dent of the fine details of the model. Unfortunately there is presently not 
enough heavy baryon data to check the formula m(£g) — m(Sg) = 0; it is not 
surprising that one gets 193 MeV instead of when one uses the hyperons 
formed with the s-quark. For the relation m(Sg) — m(Ag) « 195MeV, one 
finds 77 MeV with the ordinary hyperons but 170 MeV with the c-quark hy- 
perons. The fair agreement in the latter case is encouraging. The results in 
the former case show that the Callan-Klebanov model should not be applied 
to the strange hyperons if one goes to the true heavy meson limit. 

Additional physical quantities may be estimated with the help of the 
"potential" for the bound state problem; substituting (4.5) and (4.6) into 
(4.3) and (4.4) yields: 

Vq = -1.19(2- 0.39c + 0.26a GeV , 
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K = 0.27d + 0.12c -0.04a GeV 3 , 



(5.1) 



for the typical choice of light particle Lagrangian parameters. This, of course, 
involves the three light-heavy coupling coefficients d, c and a. The value of 
d has been discussed in the literature 21 ; bounds are obtained which agree 
with a simple estimate 16 based on pole dominance of the D — > K transition 
form factor: 



where (5 ps 1.0 and we have used 22 F D ps 250MeV. We shall adopt (5.2) for 
definiteness. In contrast we are just at the initial stage for determining the 
values of c and a; the present model may assist the direct approach from 
meson weak transition matrix elements in this regard. 

At lowest order, Vq is a measure of the energy which binds the heavy 
meson and soliton together as a heavy baryon. Allowing for a zero point 
energy Z for quantum fluctuations, we write 



The left hand side is about —780 MeV if we identify Q = b. If we first neglect 
Z and set a = 1 ( "vector meson dominance" ) , we can get an idea of the value 
of c by putting (5.1) into (5.3). This gives us c ~ 1.0, which seems quite 
reasonable in that it is expected to be of the same order of magnitude as 
d. Examining the individual terms in (5.1) we see that the pion (d) term is 
dominant but the p meson(c) contribution is 60% as large and also attractive 
while the omega (a) contribution is 40% as large and repulsive. Of course 
there can be different fits but this seems to be the typical situation. It is 
clear that the vector mesons can be expected to play a non-trivial role. 

Using the numbers given above we find the parameter k, which measures 
the strength of the R 2 part of the collective mode potential in (3.12), to be 
about 0.22 GeV 3 . With (3.12) we note that (3.11) is just the Schrodinger 



0.5 GeV. As is well known this results in a zero point energy, Z ~ 0.75 GeV. 
Such a value is not a small perturbation to the right hand side of (5.3) so 
one might expect the collective ansatz in (3.10) to be suspicious. It would be 
more reliable to include the full non-linearities of the soliton profiles as well as 
to include higher derivative interactions in the light-heavy Lagrangian. This 




(5.2) 



m(A Q ) - m N - M = V + Z. 



(5.3) 



equation for a spherical harmonic oscillator with angular 
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is a future job, but we can at least mention the parameters which will fit (5.3) 



to experiment while retaining Z — | f^) 2 - Keeping a = 1 and d = 0.53 




would result in c = 4.1 and correspondingly k, = 0.6 GeV 3 . Such a value of 
c is on the large side. An alternative fit is obtained by choosing c = 1 and 
taking d = 0.7 (the experimental bound); then we find a = — 2 (which does 
not agree with our intuitive notion of light vector meson dominance) and 
correspondingly k = 0.39. The two kinds of fits give different values for k. 
Since the orbital excitations of the spherical harmonic oscillator are equally 
spaced at intervals given by the angular frequency, we would expect the 
first orbitally excited A b to be about 0.80 GeV higher for the first fit while 
about 0.63 GeV higher for the second fit. Another possible future source 
of information about k is the Isgur-Wise function for the A& — > A c weak 
transition form factor. It was shown in ref. 9 that the non-perturbative part 
of this form factor is simply given by the overlap of the initial and final state 
spherical harmonic oscillator wave functions: 



where the strict M — > oo limit was taken for simplicity of writing. 

We would like to thank Herbert Weigel for helpful discussions. This 
work was supported in part by the U.S. Department of Energy under the 
contract number DE-FG-02-85ER40231 and also by the grant number DE- 
FG02-91ER40685. 

References 



1. C.G.Callan and I. Klebanov, Nucl. Phys. B262,365(1985); C.G. Callan, 
K. Hornbostel and I. Klebanov, Phys. Lett. B202 (1988). 

2. J. Blaizot, M. Rho and N. Scoccola, Phys. Lett B 209 ,27(1988); N. 
Scoccola, H. Nadeau, M. Nowak and M. Rho, ibid.201, 425(1988) ;D. 
Kaplan and I. Klebanov, Nucl. Phys. B335,45(1990). 

3. H.Yabu and K. Ando, Nucl. Phys. B301, 601(1988). 




r) (V ■ V) = exp \(m N 3 /K) 1/2 (V ■ V + 1)1 , 



(5.4) 



11 



4. H. Weigel, J. Schechter, N. W. Park and Ulf-G. Meissner, Phys. Rev. 
D42, 3177(1990). 

5. N. W. Park and H. Weigel, Nucl. Phys. A541,453(1992). 

6. M. Rho, D. O. Riska and N. Scoccola, Z. Phys. A341, 343(1992). 

7. E. Eichten and F. Feinberg,Phys. Rev. D 23,2724(1981); M. B. Voloshin 
and M. A. Shifman , Sov. J. Nucl. Phys. 45,292,(1987); N. Isgur and 
M. B. Wise, Phys. Lett. B232, 113 (1989); H. Georgi , Phys. Lett. B 
240, 447 (1990). 

8. E. Jenkins, A. Manohar and M. B. Wise,Caltech preprint CALT-68- 
1783-REV ( |hep-ph/9205213D ;Z. Guralnik, M. Luke and A. V. Manohar, 
Nucl. Phys. B390,474(1993); E. Jenkins and A. Manohar,Phys. Lett. 
B294,173(1992). 

9. E. Jenkins, A. Manohar and M. Wise, UCSD preprint UCSD/PTH 
92-27 (|hep-ph/9208248| ). 

10. M. Rho, Talk at Workshop on "Baryons as Skyrme Solitons", Sept. 28- 
30, 1992; Siegen, Germany (hep-ph /9210268). 

11. D.P Min, Y. Oh, B. Y. Park and M. Rho, Seoul preprint SNUTP 92- 
78 ( |hep-ph/ 92092751) ; H. K. Lee, M. A. Nowak, M. Rho and I. Zahed, 
Stony Brook preprint NTG-92-20( |hep-ph/930124l) . 

12. P. Jain, R. Johnson, Ulf-G. Meissner, N. W. Park and J. Schechter, 
Phys. Rev. D37,3252,(1988). 

13. J. Schechter, A. Subbaraman and H. Weigel, Syracuse- Tubingen preprint 
SU-4240-525( |hep-ph/9211239| ). 

14. J.D. Bjorken, in La Thuile Rencontres , 1990 , p. 583; M.Wise, Phys. 
Rev.D45, R2188(1992); F. Hussain , J. Korner and G. Thompson, Ann. 
Phys. 206,334,(1991) . 

15. M. Wise, ref. 14 above; T. M. Yan, HY. Cheng, C. Y. Cheung, G. L. 
Lin, Y. C. Lin and H. L. Yu,Phys. Rev. D46, 1148(1992); G. Burdman 
and J. F. Donoghue, Phys. Lett. B 280, 280 (1992); P. Cho, Phys. 
Lett. B 285, 145 (1992). 



12 



16. J. Schechter and A. Subbaraman, Syracuse preprint SU-4240-519 ( |hcp- 
ph/ 92Q925| ). 

17. R. Casalbuoni, A. Deandrea, N. Di Bartolomeo ,R. Gatto, F. Feruglio 
and G. Nardulli , Phys. Lett. B292, 371(1992); Phys. Lett. B294, 
371(1992);Phys. Lett. B299, 139(1993); N. Kitazawa and T. Kurimoto, 
Osaka preprint OS-GE-27-92 ( |hep-ph/9302296| ). 

18. See Appendix B of N. W. Park, J. Schechter and H. Weigel, Phys. Rev. 
D41, 2836(1990) and also Z. Guralnik et al, ref. 8 above. 

19. G. Adkins, C. Nappi and E. Witten, Nucl. Phys. B228, 552(1983). 

20. We are grateful to Herbert Weigel for supplying these. The num- 
bers given correspond to the choice of parameters 72 = 73, |<7yv</>| = 
1.9, \ h\ = 0.4, as discussed in ref. 5 above, where slightly different 
conventions from the present one were employed. Symmetry breaking 
terms were also included although their effect on the profiles is small. 

21. P. Cho and H. Georgi, Phys. Lett. B296,408(1992); J. Amundsen et. 
al., ibid, 415(1992). 

22. J. Amundsen, J. Rosner, M. Kelly, N. Horwitz and S. Stone, Enrico 
Fermi Institute preprint, EFI 92-31 ( |hep-ph/9207235l) . 



13 



